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We study the electromagnetic Dalitz decay modes of J/tp — > Pl + l~ (P = 7r°, r\ or 77'). In these 
decays, the lepton pairs are formed by internal conversion of an intermediate virtual photon with 
invariant mass 771;+;-. Study of the effective-mass spectrum of the l + l~ will shed light on the 
dynamic transition form factor Fj/^ P (q 2 ) (q 2 = m^ +; _), which characterizes the electromagnetic 
structure arising at the vertex of the transition J/tj) to pseudoscalars. We also discuss the direct 
productions of a GeV scale vector U boson in these processes J/tp — ► PU (U — > l + l~). It is 
responsible for mediating a new U(l)d interaction, as recently exploited in the context of weakly 
interacting massive particle dark matter. In this paper, we firstly use the usual pole approximation 
for the form factor to estimate the decay rate of J/tp — > Pl + l~ in the standard model. Then the 
reach of searching for the dark photon is estimated. We suggest that these Dalitz decays can be 
used to search for the light U boson in the BESIII experiment with a huge J ftp data set. 

PACS numbers: 13.20.Gd, 14.40. Pq, 13.25. Gv, 13.40.Gp, 12.38.Qk 



INTRODUCTION 

Study of the electromagnetic (EM) decays of the 
hadronic states is of interest for understanding the struc- 
ture of hadronic matter and for revealing the funda- 
mental mechanisms for the interactions of photons and 
hadrons 0,0. The EM Dalitz decays V -> Pl + l~ of light 
unflavored vector mesons (p, uj or <p) are specially inter- 
esting for probing the electromagnetic structure arising 
at the vertex of the transition V to the pseudoscalar. In 
table HI we summarize the experimental results for the 
EM Dalitz decays of the light vector mesons. The ra- 
tios of the Dalitz decays to the corresponding radiative 
decays of vector mesons are suppressed by two orders 
of magnitude, especially, for V — > Pe + e~ modes. As- 
suming point-like particles, the decay rate of this pro- 
cess vs. mi+i- can be exactly described by QED [3| in 
the standard model (SM). However, the rate is strongly 
modified by the dynamic transition form factor F\rp(q 2 ) 



_ ) , which can be estimated based on models 



of QCD BBQ 



TABLE I: The experimental results on the light vector Dalitz 
decays , V —¥ Pl + l~ (V = p, uj or <p), and ratios of the Dalitz 
decays to radiative decays of the vector mesons. These data 
are from PDG2010 \v&. 



Decay mode 


Experimental results 


r(V-tPl + l~) 
r(v->p-f) 




< 1.2 x 1(T 5 (90% C.L.) 


< 2.0 x 10^ 


uj — > 7r u e + e~ 


(7.7 ±0.6) x 10" 4 


(0.93 ±0.08) x 10~ 2 


UJ — > 7T U /i + ff 


(1.3 ±0.4) x 10" 4 


(0.16 ±0.05) x 


uj — > rie + e~ 


< 1.1 x 10~ 5 (90% C.L.) 


< 2.4 x 10" 2 


<p — > ne + e~ 


(1.15 ±0.10) x 10" 4 


(0.88 ±0.08) x 10"* 




< 9.4 x KT b (90% C.L.) 


< 0.07 x 



Experimentally, \Fvp(q 2 )\ 2 is directly accessible by 
comparing the measured invariant mass spectrum of the 



lepton pairs from the Dalitz decays with the point-like 
QED prediction. A comprehensive review of the topic 
is contained in reference [2j]. Recently, high quality data 
from NA60 experiment measured the g 2 -dependent form 
factor of the Dalitz decay us — > 7r°/i + ^ _ : 1 1[ . Using the 
usual pole approximation F(q 2 ) = 1/(1 — fa) for the 
form factor, the A~ 2 has been found to be 2.24 ± 0.06 
GeV~ 2 , which strongly deviates from the expectation of 
vector meson dominance (VMD) The form factor 

showed a relative increase close to the kinematic cut-off 
by a factor of ~ 10 For the decay of (j> ?ye + e~, 

the SND experiment has looked for the m e + e - invari- 
ant mass distribution with 213 events [l2j |. and measured 
the form factor slope A" 2 to be 3.8 ± 1.8 GeV~ 2 . Most 
recently, KLOE-2 has selected 7000 <p ~^ rie + e~ events 
with 77 — > 7r + 7r _ 7r° using a sample of 739 pb _1 on the 
4> peak 13|. A preliminary fit to the m e + e - indicates 
the possibility to reach a 5% error on the form factor 
slope fl3j ]. 

These theoretical and experimental investigations of 
the EM Dalitz decays of light vector mesons motivate 
us to study the rare charmonium decays J/tp — > Pl + l~ 
(P = tt°, r) or 77'). The measurements of the g 2 -dependent 
form factors will provide useful information on the in- 
teraction of the charmonium states with electromagnetic 
field. In particular, with more phase space, the transi- 
tion between J/yj and the pseudoscalar can be explored 
over a large region of momentum transfer, which will be 
used to test QCD prediction cleanly. The decay rates 
of the J/tp -> 77T°, 777 and 777' are (3.49±°|^) x 10~ 5 , 
(1.104±0.034) x 10~ 3 and (5.28±0.15) x 10" 3 |lo|, respec- 
tively. From a direct estimation according to the ratios of 
uj and <p Dalitz decays to the corresponding radiative de- 
cays, the expected Dalitz decay rates could reach ~ 10~ 7 
for J/tp — > 7r°e + e~ and ~ 10~ 5 for J/tp — > rje + e~ and 



respectively. 



2 



It is also interesting to search for dark photon, the 
light U boson, in J/ip —> PU and U —> l + l~, in which 
the virtual photon is replaced by the light on-shell dark 
photon U . The light U boson may couple to the SM 
charged particles with a much suppressed coupling which 
has been considered in various contexts 91 ] . We con- 
sider the new Abelian gauge group U{l)d which has a 
gauge-invariant kinetic mixing with the SM hypercharge 
U(1)y [iol-22|. After electroweak symmetry breaking, we 
have the Lagrangian: 



C = C 



SM 



e Y F 



Ad,fi Ad 



(1) 



where Csm is the SM Lagrangian, F^ v and F^ v are the 



field strength for the SM gauge boson B and U boson, 
respectively, A d is the gauge field of a massive dark U(l)d 
211 ] . The second term in Eq.([T]) is kinetic 
and e ~ 10~ 8 — 10~ 2 is generated at 



gauge group 



mixing operator, 
any scale by loops of heavy fields charged under both 
£7(l)s. In a supersymmetry theory, the kinetic mixing op- 
erator induces a mixing between the D-terms associated 
with U(l)d and U(l)y- The hypercharge D-term gets 
a vacuum expectation value from the electroweak sym- 
metry breaking and induces a weak-scale effective Fayet- 
Iliopoulos term for U(l)d- Consequently, the U(l)d sym- 
metry breaking scale is suppressed by loop factors o r by 
y/e, leading to MeV to GeV-scale U boson mass [20ll22|. 
The parameters of concern in this paper are e and rnjj . 

In the BESIII experiment, more than 1 billion J/ip (ip') 
sample will be collected next year, the first look of these 
decay modes will be accessible [23j. It will shed light on 
probing new physics beyond the standard model, such as 



possible U boson with mass less than 300 MeV range [24 - 
281 ] , which may contribute to the process by replacing vir- 
tual photon. By looking at the spectrum of the dilepton, 
with huge data set in the BESIII experiment, one may 
see possible new physics contribution which will modify 
the shape of the di-lepton spectrum and total decay rate 
of the Dalitz decay process. 

In this paper, the full angular distribution and q 2 - 
dependent rate are derived in the SM framework in Sec- 
tion n Using the pole approximation, the decay rates for 
the EM Dalitz decays of J/ip are estimated for the first 
time. An interesting consequence is that the dark photon 
can be probed in the low energy BESIII experiment. In 
Section „ we estimate the reach of the U boson search in 
the J/ip -> PU and U -> l + l~ decay. 



q 2 -DEPENDENT DECAY RATE OF J/ip -> Pl+T 

The amplitude of the Dalitz decay ip — > Pl + l~ (here- 
after, ip denotes J /ip) has the Lorentz-invariant form: 

T = 4iraf if , P e tJ ' l/f " :r p ll q I jep-^u 1 j rT V2, (2) 

where f^p is the form factor of the ip — > P transition; 
q v is the 4-momentum of the virtual photon or the total 
4-momentum of l + l~ (I — e, /a) system ; q 2 = rn 2 +l _ 
is the effective mass squared of the lepton pair; is 
the 4-momentum of the pseudoscalar meson; e p is the 
polarization 4- vector of the -0; e^ pa is the totally anti- 
symmetric unity tensor. It is straightforward to obtain 
the magnitude of the amplitude squared as below: 



irr 



167T 2 a 2 



lUp £ ){ [8(p ■ qfm 2 - 8p 2 q 2 m 2 - 2p 2 q± - 8(fc ■ p)(k 2 ■ p)q 2 + 4(p ■ q) 2 q 2 ] , 



(3) 



where m; is the lepton mass; q = k\ + fa, and p, k\ and k 2 are 4-momenta of particles P, Z + and I , respectively. 
The angular distribution of the differential decay width can be obtained as 



dT^^Pl+l-) 1 a 2 \fMl 2 )\' 



dq 2 



dQ 3 dni 



4m? 
1 



4m 



2\ 1/2 



1 - 



4m/ 

1 



[(ml -m 2 P + q 2 ) 2 - 4m 2 A g 2 ] 3/2 x 



cos 2 6»* 



(4) 



where and mp are the masses of the initial charmonium state and pseudoscalar meson; dQ^ = d^d^cosO^) is 
the solid angle of P in the rest frame of ip and dill = d^Jd(cos0J) is the solid angle of one of the lepton pair in the 
rest frame of l + l~ system (the z direction is defined as the momentum direction of l + l~ in the ip system); 9\ is the 
helicity angle of l + l~ system, which is defined as the angle between momentum direction of one of the lepton pair and 
direction of the P meson in the rest frame of l + l~ system. By integrating the solid angles in Eq. 01 one can obtain 
the g 2 -dependent differential decay width: 



dT(ip^Pl+l-) 1 a 2 \Up(9 2 )\ 2 



dq 2 



3 247rm3 



Amf 



1/2 



2mf 



[(m 



2 I 2\2 

m P -I- q ) 



4m! g 2 



3/2 



(5) 
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For the corresponding radiative decay of ip —> P7, the From Eqs. §5§ and ([6]) the q 2 -dependent differential decay 
decay width can be obtained as : width in the ip — > Pl + l~ decay normalized to the width 



1 a{mi - m 2 P f 



W Pi) = 3 - I^p(0)| 2 - (6) 



of the corresponding radiative ip — > P7 is derived 



dT{ip Pl+l~) 
dq 2 T(tp -> P7) 



a 
37 



Up(q 2 ) 



4m 2 



1/2 



\FM<?)\ 2 x [QED( 9 2 )], 



2m 2 



4m?, q 2 



3/2 



(/« 



> 2 P ) 2 



(7) 



where the normalized form factor for the ip — > P tran- 
sition is defined as F^p(q 2 ) = fjf,p(q 2 )/f^p(0), and the 
normalization is F^p(0) = 1. The form factor defines 
the electromagnetic properties of the region in which ip is 
converted into pseudoscalar. By comparing the measured 
spectrum of the lepton pairs in the Dalitz decay with 
QED calculations for point-like particles, it is possible 
to determine experimentally the transition form factor 
in the time-like region of the momentum transfer 0, Q ■ 
Namely, the form factor can modify the lepton spectrum 
as compared with that obtained for point-like particles. 

For the decays accompanied by the production of the 
electron-positron pair, we should note that the radiative 
corrections proportional to aim 2 (q 2 / mf) will be impor- 
tant. We will not discuss the high order QED corrections 
in this analysis since the data sample in the BESIII ex- 
periment is still small, and BESIII is expected to sec the 
first signal for the effect at leading order. In addition to 
that, the external conversion of the 7 from the radiative 
decay of ip —> P7 will make the analysis more compli- 
cated, however, at the BESIII the external conversion 
rate could be up to 2%, and the invariant mass of the 
m e + e - will form a narrow peak at 20-40 MeV, which will 
not really affect the slope shape of the dilepton. For the 
decays accompanied by the production of the muon pairs 
the radiative corrections and external radiation effects 
are negligibly small. 

To estimate the order of magnitude, one may use the 
Vector Dominance Model (VDM), in which the hadronic 
EM current is proportional to vector meson fields |29ll30l |. 
Hence the VDM predicts a growth of the transition form 
factors with increasing dilepton mass. The form factor 
may be parameterized in the simple pole approximation 
as 

FM1 2 ) = :— V> (8) 

where the pole mass A should be the mass of the vector 
resonance near the energy scale of the decaying parti- 
cle according to the VDM model. In ip decay the pole 



mass could be the mass of ip'. By assuming the pole 
approximation and taking A = mf, in Fig. [T] we show 
the differential decay rates for ip — > tt°1 + 1~, ijl + l~ and 
r)'l + l~, respectively. The decay rates for ip — > ir°l + l~, 
rjl + l~ and r]'l + l~ are estimated and presented in table 
HT1 To study the dependence of the decay rates on the 
value of the pole mass, we varied the pole mass. The re- 
sults for J/ip — > rjfi + are shown in Figs. [Hand [3] as an 
example. The cases for the other decay modes considered 
in this work are similar. Both the differential and total 
decay rates are not sensitive to the value of the pole mass 
if the pole mass is in the range A 2 > q 2 iax = (m.^ — m p ) 2 
in the Dalitz decay process of J/ip- The reason can be 
well understood. The dominant contribution to the de- 
cay rate comes from the region of small value of q 2 . For 
the pole mass with large value, q 2 / A 2 is small, therefore 
this term cannot give large effect. This is different from 
the case for the light vector meson Dalitz decays. 

Because the decay rates are not sensitive to the pole 
mass in the form factor, the predicted decay rates in Ta- 
ble |TT] are more reliable. Comparing these predictions 
with experimental measurement should make sense. In 
the BESIII experiment, more than 1 billion J/ip (ip 1 ) 
sample will be collected in year 2012, the first look of 
these decay modes will be accessible (23J. 



TABLE II: The estimated decay rates for ip -> tt 1 + 1~ , rjl + l~ 
and r]'l + l~ based on Eq.([7| by assuming pole approximation 
and A = m^i. The error on the decay rate is from the mea- 
sured error of ip —¥ P7 which is used as normalization. 



Decay mode 


e + e 




ip -> 7r u Pr 


(3.89^0.33) x 10-' 


(l-oili^) x 10-' 


ip -> r]l + r 


(1.21 ±0.04) x 10~ s 


(0.30 ± 0.01) x 10~ ;j 


ip -¥ r/z+r 


(5.66 ±0.16) x 10~ 5 


(1.31 ±0.04) x 10~ 5 



0.00 0.02 0.04 0.06 0.08 0.10 

m e + e - (GeV) 

(a) 




0.5 1.0 1.5 2.0 2.5 3.0 



"«,,+,,- (GeV) 

(b) 

FIG. 1: The differential decay rates for ip — > Pl + P , where the 
solid curve is for ip — > ty°1 + 1~ , the dashed curve for ip — > r)l + l~ 
and the dotted for ip — > ri'l + l~ . (a) is for the case that the 
lepton pair is e + e~, while (b) for the case . In (a), only 

part of the phase space for m e + e - is shown for demonstration 
purpose since most of the phase space is accumulated near 
small m e + e - region. 



A (GeV) 

FIG. 3: The decay rate for ip — > r/fj, + fi~ with the variation of 
the pole mass A. 



U boson as well. 

The on-shell U boson will decay to a pair of leptons 
in J/ip ->■ PU, leading to a signal of Pl+l~. The SM 
background J ftp — » P7* — > Pl + P, although large for 
this process, is not a severe problem as the kinematics 
of the signal are quite distinct. The invariant mass of 
the lepton pair is just within a single bin due to the tiny 
width of the vector U boson and can be distinguished 
from the SM background. It will be interesting to look 
for a low mass, up to GeV scale, U boson in these modes. 

From Eq. the number of background events in the 
window of Sq 2 (resolution of the q 2 ) around q 2 = mf+t- = 
mjj is about 



a. 

+ 
a. 

C- io- 



-55-1 



-a 




0.5 1.0 1.5 2.0 2.5 34) 

"V.+,/.- (GeV) 

FIG. 2: The differential decay rate for ip —t 7//i + /i~, where the 
solid curve is for the pole mass taken as A = my,/, the dashed 
curve for A = 3.0 GeV and the dotted for A = 4.0 GeV. 



REACH OF U-BOSON SEARCH IN J/ip 
U -> l+T DECAY. 



PU, 



In this section, we discuss the constraints and discov- 
ery potential for the U boson at Beijing electron-positron 
collider II (BEPCII) [Hj]. Since the U boson couples 
mainly to the SM electromagnetic current (32[. Its pro- 
duction at BEPCII is the same as that of photon, al- 
though with a much suppressed rate. Therefore, the pro- 
cess J/i/j —> P-f* — > Pl + l~ will have a chance to look for 



N B = N, 



dT(<d) -> Pl + l~) 2 

■dq 2 x BR(ip 



9mm, i 

dT(ip 



dq 2 T(iP 
Pl+l~) 



8q 2 x BR(iP ->■ P-f) 



Pi) 



V ' dq 2 T(ip -> Pi) 
N^ P (q 2 )\ 2 x [QED(q 2 )]^ 2 x BR{4, -> P 7 ), 



(9) 



where i (q^^ t ) is the upper (lower) value of the i-th 
q 2 bin; is the total number of tp decay events; and 
BR{ip — > Pj) is the branching fraction from PDG [Toj . 
The size of the bin is the window size 5q 2 which is 
obtained from the resolution functions of the m/+j- in 
Eqs. (4) and (5) in Ref [27] based on the BESIII Monte 
Carlo simulations. By replacing the photon by U bo- 
son in tp —> P7, the signal rate can be estimated to be 
BR(ip -t PU) w e 2 BR{tp -> P7), where e was defined in 
Eq. (JTJ). Thus, the expected number of signal events is 
about 



N s = N. 4 , x e 2 BR{i> -> Pj)BR(U -> l+P 



(10) 



where we assume BR(U l + l ) — 1 in this study. Com- 
bining Eqs. © and IjlOp. We estimate the expected nu- 
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merical results based on the significance 



S_ 
B 



N s 



^y/BRW ^W)BR{U -> l + n 1 , 



V^>(?)? x [QED(q2)]^2 

Therefore, with 1 billion ip events, the reach for U-boson 
searching can be e ~ 10~ 2 - 10 -3 in the ip —> PU decays. 
In Figs. [4] and [5l we show the reach of the parameter e 
by defining S/^/B = 5 for different mass of U boson. 

w 1 1 11 11 I 11 11 I 11 11 I 11 11 I 11 11 I 1 1 



y-> P eV 



10"- 



10" J r 



10' 4 r 



10"° 



0.5 



1 1.5 



2 2.5 3 

mu (GeV) 



r)'l + l~, respectively. They will be accessible in the BE- 
SIII experiment with data sample of 1 billion ip decay 
events. Especially, the g 2 -dependent differential decay 
rate can be measured by looking at the invariant mass of 
the lepton pairs. 

In the BESIII experiment, these measurements will be 
important for us to understand the interaction of vector 
charmonium states with photon, as well as to probe new 
physics beyond the standard model, we have investigated 
the signatures of a hidden U(l)d sector at the BESIII 
experiment in ip — > PU decays, and find that the BESIII 
should have an intrinsic sensitivity to the kinetic mixing 
parameter e in the range of 10~ 2 - 10~ 4 , which depends 
on the mass values of the U boson. 
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FIG. 4: Illustrative plot of the reach of vector boson at BES- 
III in the channel of ip —¥ PU (solid curve for P = 7T°; dot- 
dashed curve for P — rj and dashed curve for P = 77', respec- 
tively), where U decay into e + e~. 
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FIG. 5: The same as Fig. [4] with U decay into fi + fi 



SUMMARY 

In summary, the EM Dalitz decays of ip — > Pl + U arc 
studied in this paper. We demonstrate the differential 
decay rate as both g 2 -dependent rate and angular depen- 
dent rate explicitly. By assuming simple pole approxima- 
tion the decay rates for ip — > Pl + U arc estimated for the 
first time. The estimated Dalitz decay rates could reach 
~ 10~ 7 for ip — >• ir°l + U and ~ 10~ 5 for ip — > rjl + U and 
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